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Based on three general guiding principles, i.e., no double occupancy constraint, accurate de¬ 
scription of antiferromagnetism at half-filling, and the precise sign structure of the t-J model, a 
new ground state wave function has been constructed recently [Weng, New J. Phys. 13, 103039 
(2011)]. In this paper, we specifically study such kind of variational ground state for the one-hole- 
doped anisotropic two-leg t-J ladder using variational Monte Carlo (VMC) method. The results 
are then systematically compared with those recently obtained by density matrix renormalization 
group (DMRG) simulation. An excellent agreement is found between the VMC and DMRG results, 
including a “quantum critical point” at the anisotropy parameter a = a c ~ 0.7 (with the parameters 
t/J = 3), and the emergence of charge modulation and momentum (Fermi point) reconstruction 
at a > a c due to the quantum interference of the sign structure. In particular, the wave function 
indicates that a Landau’s quasiparticle description remains valid at a < a c but fails at a > a c due 
to the breakdown of the one-to-one correspondence of momentum and translational symmetry of 
the hole. The explicit form of the wave function provides a direct understanding on how the many- 
body strong correlation effect takes place non-perturbatively in a doped Mott insulator, which sheds 
interesting light on the two-dimensional case where the same type of wave function was proposed 
to describe the cuprate superconductor. 

PACS numbers: 71.27.+a, 74.72.-h, 02.70.Ss 


I. INTRODUCTION 


recently, whose compact form may be written as follow^® 


Ground state wave function is of great importance 
in understanding a new state of matter. Right af¬ 
ter the discovery of high-temperature superconductiv¬ 
ity in the cuprate, based on the conjecture that the 
cuprate superconductor be a doped Mott insulator, 
Anderson proposed^ a resonating-valence-bond (RVB) 
ground state, which may be simply expressed as 

I’J'rvb) =Pg|BCS), (1) 


|* G ) = ^ JVh/2 |RVB), (2) 

where |RVB) denotes a spin background (“vacuum”) al¬ 
ways remaining singly occupied. The doped holes (of to¬ 
tal number Nff) are created in pairs on such a “vacuum” 
via 


D = Y, 

ij 




( 3 ) 


where |BCS) denotes an ordinary BCS state and Pq 
the Gutzwiller projection operator enforcing the no dou¬ 
ble occupancy constraint due to strong on-site Coulomb 
repulsion. Such a ground state has been intensely 
studied 2 ® variationally since then, which will be referred 
to as the Anderson’s one-component RVB wave function. 

The no double occupancy constraint is just one of the 
most essential characterizations of Mott physics. The 
residual superexchange coupling will further cause anti¬ 
ferromagnetic (AF) correlations between the singly occu¬ 
pied spins, leading to an AF Mott insulator at half-filling 
described by the Heisenberg type model. Last but not 
least, due to strong on-site Coulomb repulsion, the orig¬ 
inal Fermi signs of the electrons will be replaced by the 
so-called phase string sign structure^®, which has been 
precisely identified in the t-J® and Hubbard 8 models for 
arbitrary doping, temperature and dimensions on a bi¬ 
partite lattice. These three constitute the basic organiz¬ 
ing principles for the strongly correlated electrons in a 
doped Mott insulator. 

Based on the above three guiding principles, a new 
class of ground state wave function has been constructed 


with the no double occupancy being automatically main¬ 
tained. Here the sign structure is implemented by a phase 
shift operator 

(4) 

i 

where 9JI ) is a statistical angle satisfying 8i(l) — 9Ji) = 
±7r with fiii defined as the down-spin number operator 
solely acting on |RVB), commuting with the c-operators 
in D. 

Here the wave function |’Fg) has a two-component 
RVB structure, with |RVB) characterizing the neutral 
spin correlation (which reduces to the true ground state 
of the Heisenberg model at half-filling) and D the Cooper 
pairing, respectively, improving the original Anderson’s 
one-component RVB state I’Frvb}- The novelty of 
Eq. § lies in that the doped holes and the spin back¬ 
ground are nonlocally entangled by the phase shift oper¬ 
ator Cli, such that each doped hole will always feel the 
influence from the background spins, and vice versa. The 
interplay will then reshape both neutral RVB and charge 
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pairing as a function of doping and result in a phase di¬ 
agram self-consistently, which provides 9 ® a systematic 
understanding/explanation of AF, superconducting, and 
pseudogap phenomena observed in the cuprate. 

In particular, if only a single hole doped into the Mott 
insulator is considered, I’I'g) in Eq. <[2j) is reduced to 

|*o)ih = £ ^(*)e-^c 4 |RVB), (5) 

i 

with tfih(i) replacing the pairing amplitude gi-j since there 
is only one hole here (created by annihilating an electron 
with 4,-spin, without loss of generality). Such a single¬ 
hole state should be contrasted with a more conventional 
Bloch type state |'I'G)ih —> |k) B L assuming iph(i)e~ zQi oc 
e’ k ri in Eq. ([H]). Here one has |k) B L —> e _ ' kl |k) BL under 
the translational transformation r; —>• r* + 1 of the hole 
coordinate in a featureless spin background. But such 
a translational symmetry of the single hole is in general 
not obeyed by Eq. ©• 

In principle, it is not a priori that the doped hole 
should carry a conserved momentum k in the Mott in¬ 
sulator, satisfying the same Bloch theorem for a doped 
hole in a semiconductor. If it does, then one says that 
the doped hole behaves like a Landau’s quasiparticle with 
well-defined charge, spin, effective mass and momentum, 
which is the basis for the Fermi liquid theory of weakly 
interacting electrons. However, the no double occupancy 
constraint in a Mott insulator means that the electrons 
are localized at each lattice site by strong interaction at 
half-filling, implying the loss of the translational symme¬ 
try of the charge. A doped hole moving on the neutral¬ 
ized spin background of thermodynamic scale does not re¬ 
store the charge translational symmetry immediately. As 
a matter of fact, in an early study of a sing le hole doped 
t-J model, it has been rigorously shown® that the hole 
acquires an irreparable many-body phase shift, i.e., the 
phase string, which demonstrates a general breakdown 
of the translational symmetry for the charge, supporting 
the argument of Anderson® that the hole doped into a 
Mott insulator does not become a well-defined quasipar¬ 
ticle due to a nontrivial scattering phase shift. 

The two-leg t-J ladder as a stack of two one¬ 
dimensional chains can serve as an ideal minimal model 
to test the novel phase string effect as well as the vari¬ 
ational wave function Eq. which can be accurately 
studied numerically by density matrix renormalization 
group (DMRG) methocP®. At half-filling, the spins are 
short-range AF correlated such that the ground state is 
gapped 13 . A single doped hole should not change the 
background spin correlation at long distance. From a 
more conventional point of view, the hole is expected to 
only carry a small distortion (spin polaron) in the spin 
background surrounding it with well-defined charge, spin, 
effective mass and a conserved momentum® 1 ' just like a 
Landau’s quasiparticle. However, the DMRG study has 
shown exotic behaviors upon doping 3 — t2 °. Specifically, 
by tuning an anisotropic paramete r a o f the two-leg lad¬ 
der, a critical value at a c is founcP^fll such that in the 


strong rung case of a < a c the doped hole indeed behaves 
like a conventional quasiparticle with a well-defined mo¬ 
mentum. However, at a > a c the momentum splits con¬ 
tinuously as a function of a — a c , accompanied by incom¬ 
mensurate char ge m odulations which violate the transla¬ 
tional symmetry 19 20 . By examining 13 19 the charge re¬ 
sponse to an inserting flux into the ring of the two-leg 
ladder, an exponential decay with the circumference of 
the ring indicates the doped charge loses its phase coher¬ 
ence or momentum conservation to become “localized” at 
a sufficiently long distance at a > a c , in sharp contrast 
to a coherent q uasip article behavior at a < a c . Further 
surprising arises 18 19 when two holes are injected into the 
gapped spin ladder, where a strong binding between the 
two holes occurs at a > a c and simultaneously the charge 
modulation disappears with restoring the translational 
symmetry. 

Microscopically, the above novel properties can all be 
attributed to the phase string sign structure hidden in the 
bipartite t-J ladder. As has b een clearly demonstrated 
in the DMRG calculation^ 13 ®®, if one artificially turns 
off the phase string sign structure in the t-J model, which 
results in the so-called a-t-J modeP®, the ordinary Bloch 
wave behavior of the doped hole is immediately recov¬ 
ered in the whole regime of a with no more critical a c . 
Simultaneously two doped holes are no longer pairecP--. 
Although how the phase string effect is responsible for 
the physics at a > a c has been qualitatively discussed in 
Ref. 1201 a microscopic and quantitative understanding of 
the DMRG results is still lacking. 

Recently, a paper by White, Scalapino, and Kivelsorl® 
has reconfirmed the existence of a c in the one-hole-doped 
two-leg ladder by DMRG under open boundary condi¬ 
tion, together with the charge modulation and the mo¬ 
mentum splitting at a > ct c found in the earlier works. 
But they gave a different physical interpretation on the 
nature of the hole state at a > a c and argued that the 
hole would still behave like a Bloch quasiparticle at mo¬ 
menta different from that at a < a c . However, many is¬ 
sues remain unanswered there, including the microscopic 
origin of the critical a c , the necessary role of the phase 
string effect that causes the charge modulation and mo¬ 
mentum splitting in contrast to the a-t-J model, the pair¬ 
ing between two doped holes, as well as the oscillation 
and exponential decay of the energy difference with the 
ladder length under a periodic boundary condition with 
inserting different fluxes, etc. In particular, how to mean¬ 
ingfully identify the Landau type quasiparticle is actually 
rather subtle in such a strongly correlated system, as to 
be clearly shown in this work. 

In this paper, we study the one-liole-doped two-leg t-J 
ladder based on the ground state I’I'g) lh [Eq. ([5])] using 
variational Monte Carlo (VMC) method. The continuous 
phase transition at a c will be naturally reproduced by 
|H'G)ih as a function of a. As a matter of fact, we find 
that a c ~ 0.7 at the t/J = 3 matches with that found 
by DMRGPMH quantitatively. Such a wave function can 
then provide a direct understanding on how a doped hole 
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moves on a Mott-insulator spin background with short- 
ranged AF correlations. 

At small a (< a c ), the ground state I'FgHii is shown 
to have a finite overlap with the simple Bloch wave state 
|k)sL at k = k 0 = (7T, 0). Conversely, if the Bloch wave 
function |k)eL is used as the variational state, the same 
k 0 is also reproduced. It means that the one-hole state 
|4 , G)ih indeed describes a coherent quasiparticle of the 
Landau’s paradigm. Namely, the quasiparticle has the 
same quantum numbers (charge, spin, and momentum) 
in both states, which differ only by the effective mass 
and thus allow for an adiabatic connection between them. 
In fact, in the limit of a —> 0, the phase factor e~ Qi 
in Eq. © can be absorbed into the hole wave function 
iph{i) such that |'Fg)iii is explicitly shown to be smoothly 
connected to the Bloch state |k)BL- 

At a > a c , an incommensurate momentum splitting 
(or Fermi point reconstruction) is exhibited in |’I'G)ih> 
accompanying a charge density modulation, also in ex¬ 
cellent agreement with the DMRG result. However, if 
one still uses the Bloch state |k)eL as a variational state, 
the momentum carried by the hole is found to always re¬ 
main commensurate at ko = (0,0), such that there is no 
overlap between I'FcHh and |k 0 )BL at a > a c . Conse¬ 
quently the doped hole can no longer be described as a 
Landau-type quasiparticle due to the breakdown of the 
one-to-one correspondence between a bare electron and 
a quasiparticle with the same momentum. 

Here the incommensurate momenta in |\&G)ih as well 
as the charge density modulation can be clearly related 
to the intrinsic quantum interference pattern due to the 
phase string effect. The results are in sharp contrast 
to the scenario in which the charge modulation is sim¬ 
ply interpretecf 21 ! as a standing wave of two opposite- 
propagating Bloch waves mixed only by the reflection 
at open boundaries of the two-leg ladder. In other words, 
the present charge modulation is related to the absence 
of the translational symmetry in the bulk of the ladder 
where the phase string effect is unscreened at a > a c 
because of the spatial sep aratio n of the hole and its spin 
partner with increasing a 19 20 . On the other hand, in 
the a-t-J model where the phase string sign structure is 
absent, we show that the single-hole ground state indeed 
reduces to the Bloch-like one |k 0 )BL with k 0 = (0,0), 
in which the Landau’s quasiparticle description works in 
the whole regime of 0 < a < 1, again in good agreement 
with the DMRG result. 

The bottomline is that the ground state |4 f c)ih can 
well capture the essential properties of the one-hole- 
doped two-leg t-J ladder found in the DMRG simula¬ 
tion. Given the gapped spin vacuum, these novel proper¬ 
ties are solely associated with the phase string induced by 
the doped hole, whose effect cannot be reduced to simply 
renormalizing the effective mass of a Landau’s quasipar¬ 
ticle at a > a c . Here the many-body phase string factor 
Cli in Eq. ([H]) prohibits a perturbative approach starting 
from a Bloch state, leading to the intrinsic translational 
symmetry breaking of the charge. In the end of the paper, 


we shall also briefly discuss how the phase string effect 
further renders the hole self-localized through an ulti¬ 
mate translational symmetry breaking, which involves a 
many-body correction to the hole wave function tph(i)- 
The rest of this paper is organized as follows. In Sec.|TTJ 
we introduce the model and further discuss the Bloch-like 
and non-Bloch-like single-hole variational wave functions. 
The procedure of the VMC for determining the varia¬ 
tional parameters of the wave functions is then outlined, 
with the details presented in Appendix |Aj In Sec. m 
we identify a second order phase transition at a c . The 
changes of physical properties, such as charge density 
modulation, momentum distribution and quasiparticle 
weight, are investigated by VMC and compared with the 
DMRG results. Thereafter, the physical nature of the 
variational wave function is further examined. As a com¬ 
parison, we then show both analytically and numerically 
that the Bloch-like wave function well captures the es¬ 
sential properties of the a-t-J model, of which the sign 
structure is merely the Marshall sigrP 2 rather than the 
phase string (see Appendix |B]) . Finally Sec.[V]is devoted 
to the conclusion and discussion. 


II. VARIATIONAL APPROACH 

A. The model 

In this work, we focus on a n anisot ropic t-J model on 
a two-leg square lattice laddei^GSESHm 

Hu = E ( H h + H i) > (6) 

<y> 

H ij = ocijt Y c\ a c jt7 + h.c., (7) 

a 

= a.ijJ (s 4 • Sj — —, (8) 

where a, j = a (aij = 1) if the nearest-neighbor bond (ij) 
is parallel (perpendicular) to the chain direction. The di¬ 
mensionless parameter a controls the anisotropy of the 
model. We fix t/J = 3 in this paper for simplicity. Here 
Si and n t are the electron spin and number operators at 
site i, respectively, with the no double occupancy con¬ 
straint rii < 1 always enforced on the Hilbert space. 

We shall also discuss the so-called a-t-J model 13 -E^H 20 ] 
on a two-leg ladder lattice, which differs from the t-J 
model only by replacing the hopping term if *• in Eq. (|7|) 
with 

Hij 1 = -Oiijt Y ac L c j° + hx -< (9) 

<7 

where a spin-dependent sign a = ±1 is added to each 
step of hole hopping. A comparative study of these two 
models will reveal the fundamental physics hidden in the 
Mott physics. 
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B. Ground state at half-filling 


At half-filling, the two-leg t-J ladder reduces to a pure 
Heisenberg spin ladder, whose ground state |RVB) is AF 
short-range-correlated, se parat ed by a finite energy gap 
from the first excited stat d 13 * 1 ^ . As a good starting point, 
Liang-Doucot-Anderson type®^ bosonic RVB variational 
wave function will provide an excellent description: 


|RVB) = w v\v). 


( 10 ) 


Here \v) = £ w (il(y)eo cU-cj^jO) is a sin¬ 

glet pairing valence bond (VB) state specified by the 
dimer covering configuration r 3 '. The antisymmetric 
Levi-Civita symbol e CTi>(Tj . ensures the singlet paring be¬ 
tween spins on sites i and j. The amplitude of each VB 
state |u) is factorized by w v = II(»j)ei; ; w here is a 
non-negative function depending on sites i and j belong¬ 
ing to opposite sublattices, respectively. Such a factor¬ 
ization tremendously decreases the number of variational 
parameters, with hij ’s chosen such that (RVBjRVB) = 1. 


Moreover, the variational wave function Eq. (101 satis¬ 
fies the exact Marshall sign rule for bipartite H< 
rnodel^. 


eisenberg 


By using the wave function |RVB) in Eq. (101, we 
perform the VMC calculation^® of the superexchange 
Heisenberg model Eq. ([8]) on a lattice with size N = 
N x x N y = 40 x 2 under an open boundary condition. 
The optimized superexchange energies Ej calculated by 
VMC are in an excellent agreement with the DMRG en¬ 
ergies for different a as shown in Fig. 0(a)- 


C. Single-hole-doped variational state 



FIG. 1. (color online). The superexchange energy Ej shows 
no singularity as a function of the anisotropic parameter a. 
(a) The VMC calculation for the state |RVB) in Eq. (101 (red 
crosses) and the DMRG result (black dots) in the pure two-leg 
spin ladder; (b) The superexchange energy of the single-hole- 
doped variational state Eq. © and the DMRG result. The 
two-leg ladder has size N = N x x N y = 40 x 2 under the open 
boundary condition. 


Once the half-filling ground state |RVB) is accurately 
known, the ground state for the single hole case can 
be generically constructed by removing an electron (say, 
with a 4- spin) from the vacuum state |RVB) as follows 

|*G>ih = l>h(i)c4|HVB> . (11) 

i 

The wave function &h{i) is a many-body operator which 
involves the hole coordinate i and at the same time can 
be spin-dependent. Generally the latter accounts for the 
spin background response to the creation of a bare hole at 
site i, cal RVB). It is usually called the “spin-polaron” 
efFectl. 14 17 . A quasiparticle description is valid if such 
an effect remains featureless, meaning that it only renor¬ 
malizes the effective mass and the wave function spectral 
weight without changing the momentum. In particular, 
in the present two-leg ladder case, spins are gapped at 
half-filling over the whole range of 0 < a < oo such that 
the correction of the spin-polaron effect to the renormal¬ 
ization is expected to be weak. Thus, in the present work 


we shall always neglect a featureless spin-polaronic cor¬ 
rection to $;,(*), although it may be still important to 
improve the variational energy. 


1. A Bloch-like state 

Then, if the whole spin-polaron correction to the bare 
hole “Wannier basis” CjjJRVB) is neglected, <1V(*) will 
reduce to a single-particle Bloch wave function 



by assuming a translational symmetry for the hole (which 
is not a priori in a many-body system, see below). Cor¬ 
respondingly |^G)ih = |k)eL is uniquely specified by the 
momentum k of the hole without involving any other 
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variational parameters: 


k >BL = \/^E e,;kriC 4|RVB>. 


(13) 


2. A non-Bloch-like state 

Nevertheless, even if the longitudinal (amplitude) spin- 
polaronic effect is negligible in a spin gapped background, 
a transverse or many-body phase shift of the spin back¬ 
ground in response to the creation of the bare hole may 
still play a crucial role in the present strongly correlated 
system. Specifically, one may construct a new variational 
wave function as given in Eq. ([5|: 

®h(i) = Vh{i)e~ lhi , (14) 

where the phase factor , defined in Eq. Q, is 

a nonlocal operator depending on the spin configura¬ 
tion in the vacuum. [The normalization (’I'gI’I'g) = 
1 implies the normalization of the hole wave function 
J2i Wh{i)\ 2 = 2. ] Note that the new “Wannier basis” 
cqjRVB) = e -lf2i cq.|RVB) still remains invariant under 
the whole hole-spin translational operation. But <fih(i), 
determined variationally as a single-hole wave function, 
is no longer necessarily Bloch-wave-like as in Eq. ©■ 

It is important to point out that, in contrast to a con¬ 
ventional weakly-interacting system, a Bloch wave con¬ 
struction is not automatically valid for a strongly corre¬ 
lated many-body system. Here one actually deals with 
a doped hole moving in a quantum spin vacuum rather 
than an inertia translationally invariant vacuum, say, in 
a semiconductor. In the former, the translational sym¬ 
metry for the charge is not upheld generally for a relative 
motion with regard to the charge neutral (Mott insulator 
) spin background. (Note that this does not contradict to 
the translational symmetry of the total system composed 
of the charge and spins as a whole.) As a matter of fact, 
it was rigorously shown 5 that a hole transverses along 
a closed path c in a doped Mott insulator, described by 
the t-J model on a bipartite lattice, will always pick up 
a nontrivial phase string factor (—l) N ±( c \ where IV^(c) 
denotes the total number of down spins exchanged with 
the hole along the path. Clearly (—l)Ni.( c ) represents a 
non-integrable (path-dependent or Berry-like) phase fac¬ 
tor associated with the motion of the doped hole, which 
generally breaks the translational symmetry. 

In the present single-hole-doped two-leg ladder, one 
may expect that the phase string effect be strongly re¬ 
duced over a long-wavelength scale due to the presence 
of an energy gap in the spin background, in contrast to 
a gapless case. However, due to the singular and non¬ 
local nature of the phase string picked up by the doped 
hole, it is still very crucial to carefully treat such an effect 
in an energetic variational procedure, which involves the 
nearest-neighbor hopping and superexchange processes 


where the quantum interference of the phase strings plays 
a critical role. 

The phase-string operator Cli in Eq. ([ 4 ]) can produce a 
phase shift ±7r each time when the hole and a down-spin 
exchange positions during the hopping. In this way, the 
above-mentioned singular phase string (—1)^ ^ gets ac¬ 
curately encoded by e~ l ^ li in Eq. Consequently (fih(i) 
becomes a much smoother wave function which can be 
then determined variationally. In this sense, the phase¬ 
string factor e~ lQi regulates the singular phase string 
effect in the t-J model and transforms the model into 
a perturbative-treatable formalism (its version at arbi¬ 
trary doping [cf. Eq. ([ 2 ])] has been previously obtained 
in Ref. [S]). Here it is instructive to point out that the 
topologi cal ph ase string factor e _if2i enforces the mutual 
statistic. j 38 l 39 l between the doped hole and the down-spins 
as indicated by the above sign structure. It plays the 
same role as the statistical phase factor Y\ i<: j{zi — Zj ) 3 in 
the Laughlin wave function for v = 1/3 fractional quan¬ 
tum Hall systerrPni, which ensures the anyonic statistics 
of the same species r ather than the two different species 
as in the present cas d 38 * 39 ( In both cases, a traditional 
perturbative analysis is applicable only after explicitly 
identifying the topological/statistical phase factor. 

The statistical angle 9i(l) may have different choices. A 
natural and symmetric choice of Qi(l) in two-dimensions 
is 9i(l) = 9®[l) = lining — zi), where Zi = Xi + iyi is the 
complex coordinate. In one-dimension, it further reduces 
to 9i(l) = ±7r 9{i — l) with 9(i — l) = 1 or 0 at i > l 
or i < l , respectively. In the present anisotropic two- 
leg ladder, one may introduce a variational parameter A 
(0 < A < 00 ): 

tan 9i(l) = A tand 9 (l) , (15) 

where e l9i ^ and e l8 i9) are located in the same quadrant 
of the complex plane. In the following we shall see that 
A —i 0 in the strong rung regime a < a c , while A —> 00 
in the decouple chain limit a —> 00 . 


D. Variational procedure 


Based on the variational wave functions given in 
Eqs. ([5]) and (131, one can decide the ground state by 
optimizing the total energy via a VMC procedure out¬ 
lined as follows. 

(1) Firstly, the half-filling ground state |RVB) in 
Eq. ( |To] ) is optimized as discussed in Sec. IIB [cf. Fig. [T] 
(a)]. Upon doping one hole into the two-leg spin lad¬ 
der, the variational parameters fty-’s should remain un¬ 
changed in the thermodynamic limit. 


(2) Based the Bloch-like wave function Eq. (131, one 
finds that the hopping energy is given by 


E t = {H t )bl = 2fsL cos + tgL cos k v 


(16) 


with tgL = t( 1 +4(RVB|S i ■ S ,-|RVB))/4 (ij are the near¬ 
est neighbors). From Eq. (16), one sees that the only 
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variational parameter is the momentum k which mini¬ 
mizes the hopping energy at kn = (w, 0) if fS T > 0 and 
k 0 = (0,0) if f| L < 0 (with t^ L < 0). 

(3) On the other hand, based on the non-Bloch-like 
wave function Eq. ([5]), 


E t = (H t ) = - 'Y^i ij V* h {j)‘P h {i) + h.c., (17) 

<*i> 


in which the hopping matrix element is given by 


tij — Otij ^, 


(RVB|cT ; c J - <r e i ^- <A ‘ct r c 4 |RVB) , (18) 


which can be directly computed in the variational calcu¬ 
lation (see Appendix A3). One then determines iph{i) by 
diagonalizing Eq. (|17p under a given A. Finally, the total 


energy E to t = E t + Ej is minimized by optimizing A. 

We shall carry out the above-sketched variational 
procedure by using the loop update Monte Carlo 
algorithirP^. The detail formulas used in computing the 
above superexchange and hopping energies are given in 
Appendix |A2| and |A3[ respectively. 


III. GROUND STATE PROPERTIES OF THE 
TWO-LEG LADDER DOPED BY ONE HOLE 



a 


FIG. 2. (color online). Hopping energies Et(a) for the non- 
Bloch-like wave function Eq. ([S]), the DMRG ground stat^^, 
and the Bloch-like wave function Eq. ( |13| ), respectively. The 
inset shows the second order derivative of the DMRG energy, 
indicating a second order transition takes place at a c ss 0.7. 
A similar critical point of a c ~ 0.7 is also identified for the 
non-Bloch-like wave function Eq. ©■ By contrast, there is a 
level crossing at a® L ~ 0.6 for the Bloch-like wave function 
Eq. (13 1 with the momentum shifted from ko = (tt, 0) to 
(0, 0). Note that the above VMC and DMRG calculations are 
carried out at a finite ladder size N = 40 x 2 under the open 
boundary condition. 


Based on the VMC calculation outlined in the previ¬ 
ous section, we present the ground state properties of 
I^G>ih m Eq. 0 below, in comparison with the DMRG 
simulations as well as the conventional Bloch state jk)uL 
satisfying the translational symmetry. 


A. “Quantum critical point” at a = a c 

The ground state energy of the two-leg t-J ladder, E tot , 
is composed of the hopping energy E t and the superex¬ 
change energy Ej. As the starting point at half-filling, 
|RVB) gives rise to an excellent energy Ej in comparison 
with the DMRG result as shown in Fig. [I] (a). Ej remains 
a smooth function of a for |^c)ih upon one-hole-doping, 
which is shown in Fig. |T] (b) together with the DMRG 
data. 

On the other hand, according to the DMRG 
calculatiorP-^, the hopping energy E t of the single hole 
shows a “quantum critical point” at a c ~ 0.7 as indi¬ 
cated by the second derivative over a (see Fig. [2] and the 
inset). By comparison, the corresponding kinetic energy 
of |R/g) ih is also shown in Fig. |2j which exhibits a singu¬ 
larity at a c « 0.7 (see below) indicated by the red arrow, 
very close to that of the DMRGP^. 

To look more closely, in Fig. i Etot of |^G)ih as a 
function of the variational parameter A is presented at 
two typical values of a (a = 0.3 < a c and a = 0.8 > a c ). 
One finds that the energy minimum takes place at A ~ 0 
if a < a c , and A ^ 0 if a > a c . The inset of Fig.[3]further 



FIG. 3. (color online). The critical transition at a c deter¬ 
mined by the non-Bloch-like wave function Eq. as shown 
in Fig. [5] When a < a c (a > a c ), the minimal total energy 
Etot corresponds to A —► 0 (A 0). The inset shows the “or¬ 
der parameter” A as a function of a [E to t(A = 0) is rescaled 
to zero] under a finite ladder size N = 40 x 2. 


shows A vs. a. The systematic change of E to t(\) with 
respect to a thus resembles the Ginzburg-Landau theory 
of second order phase transition, with a playing the role 
of the “temperature” and A the “order parameter”. Since 
a finite size calculation is involved here [Fig. [3], one needs 
to examine more carefully the distinct behaviors on the 
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two sides of a c in the following. 

In contrast to the continuous transition at a c found in 
the DMRG and the ground state |d/G)ih, the kinetic en¬ 
ergy of the Bloch state |k)sL [defined in Eq. (13)] shows 
instead an abrupt change (level crossing) from the mo¬ 
mentum k 0 = (n, 0) to k 0 = (0,0) at a)? L ss 0.6 with the 
increase of a, which is presented in Fig. @by the dashed 
curves. Such a first order transition is simply due to the 
sign change of the effective hopping parameter tg L with 
the decrease of (RVBjSi • S.j+^jRVB) in Eq. (16). 




k/jt 


FIG. 4. (color online). Typical hole momentum distribution 
1 — 7ik for a < etc and a > a c : (a) the VMC calculation; 
(b) the DMRG calculation. Here each peak in (a) and (b) 
denotes a characteristic momentum, which is split from the 
commensurate momentum ko = (w, 0) at a < a c to k x = ivdcK 
( k y = 0) at a > a c . 


B. Bloch-wave behavior at a < a c 

Let us examine the nature of physics on the two sides 
of a c in detail. It has been found by DMRGP^ that the 
hole momentum distribution 1 — nk (n-k = ]Cer( c ka- c k<7)) 
is peaked at momentum k 0 = (7r,0) at a < a c , which 
is then split into two peaks at a > a c . Very similar 



FIG. 5. (color online). Quasiparticle weight at a = 0.5 < 
a c and a = 0.8 > a c , respectively. The peak position indi¬ 
cates that the characteristic momentum is shifted from the 
commensurate ko = (7r, 0) at a < a c to k x = n ± k ( k y = 0) 
at a > a c . The inset shows that Z to t = X)k -^ k * s sm °°th and 
nonzero across a c as a function of a (see the text), which is 
consistent with the DMRG result 2 ®. 


properties are found for |H/G)ih in the VMC calculation 
[see Fig. [4] (a)], which are in good agreement with the 
DMRG results [Fig. [4] (b)]. 

It indicates that at least at small a (< a c ), |'I'G)ih an d 
|ko)BL may describe the same quasiparticle state. Note 
that the main distinction between the variational wave 
function Eq. © and the Bloch wave function Eq. ( fTTT[ ) 

lies in the many-body phase factor appearing in 

the former. To compare these two wave functions, we 
study the wave function overlap defined by 


Ok = BL(k|^ r G)lh 


(19) 


= \-y 

V N 


— ik-r* 


^(f)(RVB|e- in ’n 4 |RVB). 


Correspondingly the “quasiparticle spectral weight” is 
defined by Z k = |ak| 2 , which measures the probability 
of finding the bare h ole st ate |k)sL m the ground state 
I^G)ih (see Appendix |A5| ). Then the ground state of the 
variational wave function Eq. (14) may be reexpressed as 
follows 


|^G)lh = «k|k)BL H-) (20) 

where the second term • • • on the right-hand-side (rhs) 
refers to the non-Bloch-like part that is orthogonal to the 
bare hole state in the first term. 

We have shown that using the first term alone, i.e., 
the bare hole state which is a Bloch-like state, in the 
variational procedure, will result in a commensurate mo¬ 
mentum at ko = (7r, 0) at small a (< a c ). The two wave 
functions will thus have a finite overlap Z k at k = k 0 (cf. 
Fig. [5]), implying that |4/G)ih has the same momentum 
























0.08 


(a)VMC —■—a=0.5 




FIG. 6. (color online). The hole density n = \ifih(i )\ 2 /2 is 
smooth for a < a c , and oscillating for a > a c under an open 
boundary condition. Note that ri( l = n ( ' x y j is independent of 

y and n h (x) = T, y n (x, y )- 


1.2 - (a) VMC 




FIG. 7. (color online). Two characteristic momenta Q o and 
2 k turn out to coincide with each other: Qq is the charge 
modulation momentum; 2k is defined as the double-peak dis¬ 
tance in rik or Z k. They vanish for a < a c and approach n 
for a a c - 


k 0 . Indeed, our VMC calculation shows that 

(RVB|c^e- lf2 -c 4 |RVB) oc e ±lk ° r ‘ (21) 


C. Charge modulation as the fingerprint of 
translational symmetry breaking at a > a c 


[which may be understood analytically as A —> 0 in 
Eq. ^] and <Ph(i) ex constant to result in Z^ 0 ^ 0. 
In this regime, the distinction between |'Fg) and |ko)BL> 
i.e., the second term on the rhs of Eq. (201, is mainly re¬ 
sponsible for the effective mass and kinetic energy renor¬ 
malization without changing the momentum k 0 . 

Thus, at a < a c , even though the ground state energy 
may get further improved by the phase string factor 
as shown in Fig. [2j the Bloch wave description Eq. (13) 
still remains qualitatively valid with a correct momentum 
k(j. It is consistent with the general Landau’s paradigm 
that the quasiparticle wave function has a finite over¬ 
lap with the Bloch wave function of a bare hole, sharing 
the same quantum numbers including the one-to-one cor¬ 
respondence of the momentum. As emphasized before, 
the featureless spin-polaron effect, which is neglected in 
I’I'g), nray improve further the kinetic energy, but is not 
expected to change the above one-to-one correspondence 
of the momentum at ko- 


The translational symmetry is underlying the above- 
discussed Bloch-wave-like description of the doped hole 
at a < a c . But such a symmetry will be found broken at 
a > a c in the non-Bloch-like wave function Eq. |5]). 

By VMC one finds that the hole density n * is smooth 
at a < a c , but becomes oscillating at a > a c as illus¬ 
trated in Fig. [b] (a). Here the hole density nh can be 
related to the variational hole wave function as fol¬ 

lows (see Appendix |A 4[ ) 

n i = 1 “ X^ C !<7 C ^> = \ \<Ph(i)\ 2 - (22) 

G 

The DMRG results are shown in Fig. [6] (b) for the same 
parameters, and one finds that the VMC and DMRG are 
in a qualitative agreement. 

The corresponding charge modulation wavevector Q o 
is shown in Fig. [7] (a), which vanishes at a —> a c and 
approaches n in the large a limit. Here Q o well matches 
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with the momentum shift 2 k between the two peaks in 
Fig. m where the original momentum peak at k = ko 
is split into double peaks at k± = k 0 ± (k, 0 ) at a > 
a c . These results are once again well consistent with the 
DMRG resultiPSI presented in Fig. @(b). 

The charge modulation wavevector Q o can be therefore 
used to quantify the qualitative change (“phase transi¬ 
tion”) of the ground state |'I , G)ih at ct c , which may be 
more physical than the variational parameter A shown in 
Fig- 0 Here the spatial oscillation of the charge density 
can be traced back to the flux structure of tjj defined in 


the step-(3) of the variational procedure in Sec. IID [cf. 
Eq. @]. At a > a c , the mean-field solution of ifihJT) be¬ 
comes oscillating with breaking translational symmetry 
to imply the unscreened phase string effect (see below). 


D. Breakdown of Landau-type quasiparticle 
description at a > a c 

As already seen previously, the overlap between |’Fg)iii 
and |ko)sL disappears, i.e., Z^ a = 0 at a > a c because 
the momentum in the former is split into incommensurate 
peaks k± which no longer coincide with the commensu¬ 
rate k 0 at (0,0). According to Z k shown in Fig. [5| the 
momentum is shifted to the incommensurate positions 
k ± instead. Here the distinction between |'F g ) 1 i 1 and a 
Bloch state with momenta shifted to k± is crucial. In 
order to get the correct momentum k-t, one cannot sim¬ 
ply start with the bare hole or the Bloch state |k)sL in 
the first term of Eq. (20). As noted above, one would 
always find a commensurate ko if the Bloch state jk)eL 
is to be used alone variationally or as a self-consistent 
mean-field state. Rather one has to utilize the full form 
of I’FgHii including the non-Bloch-term denoted by • • • 
on the rhs of Eq. (20). Therefore it is no longer pos¬ 
sible to “adiabatically connect” |'I / G)ih with |k±)BL at 
a > a c , because the latter cannot get the energy (in¬ 


volving nearest-neighbor hopping process) nor momen¬ 
tum (long-wavelength physics) right as a stable rnean- 
field/variational state. This clearly signals the break¬ 
down of Landau’s one-to-one correspondence assumption 
of the momentum for the quasiparticle. 

It is instructive to further examine how the Landau- 
type quasiparticle picture breaks down even though 
Ztot / 0 as shown in the inset of Fig. [5j which does 
not show any singularity at a = a c , consistent with the 
DMRGP. Here Ztot is defined by Ztot = Ek Zk. — 
E i ](RVB|c 4 |tf G )] 2 , which measures the probability of 
the true ground state remaining in a bare hole state. It 
can be further expressed as 

Ztot = E I^W| 2 |(RVB|c| i e- i ^c 4 |RVB)| 2 . (23) 


In evaluating at in Eq. (19) or Ztot in Eq. (23), e 


n—iQi 


is averaged over the half-filling state |RVB), which gives 
rise to a trivial numerical oscillator similar to Eq. (21) 
even at a > a c to result in a finite Z to t- 


On the other hand, the incommensurate splitting of 
k-t is decided by the wave function iph(i) as the solution 
of Eq. (171 in the variational procedure. Note that in 
Eq. (171, the effective hopping integral L, in E q. (|18| ) 
may be expressed analytically via Eqs. ([5j and (|14|) as 
iij = (RVB|tjj|RVB), with 


2 — J.0 TT 

L ij ~ L ij- rL iji 


(24) 


where ~ a ijt (j — |(RVB|Sj ■ S ? |RVB)) is accompa¬ 
nied by a phase factor 


Hij = e l ( Ai i ^*2. 


(25) 


Here the phase-string factor comes into the crucial 

play: its phase difference during the nearest neighbor 
hopping gives rise to a nontrivial flux per plaquette in 
Hij via the gauge link variable^ 


A h - < = E ^ 0 i(O)«4 


(26) 


with 7iq = cj cij. In the present variational approach, 
one can numerically determine the flux associated with 
tij with the solution <Ph{i) exhibiting the charge modu¬ 
lation shown in the last subsection. 

Upon a careful examination, one finds that the effec¬ 
tive flux will be sensitive to the spins near the hole, as¬ 
sociated with the hole creation in e,|RVB), as the rest 
of the spins are in the short-range RVB paired vacuum 
whose contribution to effectively diminishes away 

from the hole site. At a <C 1, the RVB pairs are mostly 
rung-paired such that the spin partner of the doped hole 
is sitting at the same rung of the hole. In this limit, one 
has A —> 0 and the effective flux vanishes to result in a 
translational invariant state. At a larger a, the separa¬ 
tion between the hole and its spin partner gets enlarged 
so that becomes nontrivial with A / 0 (i.e., the 

phase string becomes unscreened—*^) to result in the 
new phase at a > a c . 

Therefore, the single-hole ground state is indeed cor¬ 
rectly described by the variational wave function Eq. ([5]) 
rather than the Bloch-like one Eq. (13) at a > a c , 
where the phase-string factor e~ lQi is crucial in regulat¬ 
ing the singular short-range hopping process to optimize 
the ground state energy. In particular, the incommensu¬ 
rate splitting in momentum k± is decided by <Ph (*) with 
g-iO; paying an indispensable role in the variational so¬ 
lution. Furthermore, due to the flux effect associated 
with the effective hopping integral tij, the single-hole’s 
translational symmetry is generally broken. 


E. Beyond the simple variational theory: 
Localization 


At a > a c , with the breakdown of the one-to-one 
correspondence, the incommensurate momenta k±’s are 
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FIG. 8. (color online). Absence of the phase string sign structure: VMC and DMRG results for the a-t-J model, (a) Hopping 
and superexchange energies as a function of a. (b) Hole density rfl shows no charge modulations at a = 1. (c) Momentum 
distribution 1 — has a peak at k = (0, 0) at a = 1. 


no longer “protected” as in the Landau’s quasiparticle 
description, which may be subject to fluctuations be¬ 
yond the variational/mean-field theory. To go beyond 
the variational or mean-field approach, one may gen¬ 
erally consider <Ph(i) as a many-body wave function in 
Eq- 0- Then the hopping energy ( H t ) may be rewritten 
as (H t ) = (RVBI^IRVB) with 

H h = ~^2 iij<Phti) l Ph(i) + h.c. (27) 

(ij) 


Compared to the previous variational scheme in deter¬ 
mining the single-hole wave function iph(i) i n Eq. ( |17| , 
now the wave function tph(i) becomes a many-body one 
that is directly subject to the gauge flux [Eq. ( [26] )] ap¬ 
pearing in t„j :l before the average over |RVB) is taken. 
The flux enclosed within a plaquette as contributed by 
Eq. (26) is given bjll 


- C-) = n Yl ( n 4 -• ( 28 ) 

□ lea 


Besides an average flux 0o, one may estimate the fluc¬ 
tuation S(f >□ =0o—0o by (RVB|(50 d) 2 |RVB). Since 
this is a quasi-one-dimensional system, such short-range 
plaquette flux fluctuations at A ^ 0 are generally ex¬ 
pected to cause the Anderson-like localization of the 
doped hol^2l. In other words, any unscreened scat¬ 
tering between the two characteristic momenta k± at 
a > a c , as caused by 5<f>n, will inevitably lead to the 
self-localization of the charge. This is consistent with 


the DMRG calculatioi i — 112 1 for the two-leg ladder un¬ 
der t he per iodic boundary condition. There it has been 
showrpUnD that the energy difference caused by the in¬ 
serting flux into the ribbon shows an oscillation and expo¬ 
nential decay with the ladder length. Such an insensitiv¬ 
ity of the doped charge to the inserting flux indicates that 
the hole becomes phase incoherence in a sufficient large 
ladder or self-localized as far as the external U(l) gauge 
field is concerned. The VMC based on Eqs. ( [27] ) and (28) 
has indeed confirmecP^ the charge localization a > a c in 
agreement with the DMRG. The self-localization of the 
hole and the detailed behavior of Z ^ near k-t in the large- 
ladder-length limit will be further discussed elsewhere. 


IV. a-t-J MODEL: ABSENCE OF THE PHASE 
STRING 

So far we have focused on the t-J model. The phase 
string sign structure has been considered to be the most 
essential factor, which gives rise to the phase string oper¬ 
ator in the ground state Eq. Now we consider 

the case where such a sign structure can be precisely 
removed^, which results in a modified local Hamilto¬ 
nian known as the cr-t-J model with a distinct hopping 
term Eq. Then the distinction between the t-J and 
a-t-J models can directly tell the singular role of the 
phase string effect, which has been clearly demonstrated 
by DMRG simulationd 1 ®^®! as mentioned before. 

As shown in Appendix [B] the sign structure of the 
a-t-J model can be rigorously identified as the Marshall 
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sigrP2] means that after a Marshall-sign transforma¬ 
tion, the a-t-J model (in arbitrary dimension and hole 
concentration) can be transformed to a model with triv¬ 
ial sign structure as defined in Ref. 1441 According to 
the Perron-Frobenius theorem, the ground state of this 
model has non-negative coefficients in the basis satisfy¬ 
ing the Marshall sign rule^l In particular, the Bloch-like 
wave function Eq. (13) with ko = (0, 0) satisfies this sign 
structure. 

It means that the one-hole-doped ground state of the 
a-t-J ladder should be well described by the Bloch-like 
wave function Eq. (13), in which the translational sym¬ 
metry is expected to be generally maintained (the fea¬ 
tureless spin-polaron effect is still negligible due to the 
spin gap in the background). Fig. [8] clearly illustrates the 
overall agreement of the Bloch state |k 0 )BL at k 0 = (0, 0) 
with the DMRG result for the one-hole-doped a-t-J lad¬ 
der. 


V. CONCLUSION 

A single hole injected into a two-leg spin ladder has 
manifested a series of novel prop erties as recently re¬ 
vealed by DMRG simulation :) 1 3 b In this work, we 

have studied such a system based on a variational wave 
function |'kc)ih hr Eq- (5| using VMC method. An ex¬ 
cellent agreement with the DMRG results have been ob¬ 
tained, which suggests that the trial ground state Eq. |5| 
has captured the most essential features of such a doped 
Mott insulator. 

The foremost important message delivered in this work 
is that the phase string sign structure plays a critical role 
in a doped Mott insulator. Indeed, by artificially switch¬ 
ing off the phase string sign structure in the t-J model to 
result in the a-t-J model, both DMRG and VMC calcu¬ 
lations have shown that the exotic properties exhibited in 
the former model are totally replaced by a conventional 
Bloch-wave behavior of the doped hole similar to that in 
a translationally invariant semiconductor. 

In essence, the nontrivial phase string effect implies 
the translational symmetry breaking in a doped Mott 
insulator. Both DMRG and VMC have shown that such 
an effect is responsible for the emergent critical point a c 
in the anisotropic two-leg t-J ladder of the single hole 
case. The translationally invariant Bloch state of the 
doped hole only survives in the strong rung regime of 
a < a c , where the phase string gets “screened” with 
A —► 0 due to a tight binding of the hole with its spin 


partner moving in the spin gapped vacuum? 19 -®. The 
exotic phenomenon arises at a > a c where the phase 
string starts to become “unscreened” with A ^ 0 as the 
separation between the h ole a nd its spin partner gets 
enlarged with increasing a ! 19 l 20 l . 

The fingerprint of the unscreened phase string is char¬ 
acterized by the emergent charge density modulation at 
a > a c . Based on |4'G)ih, one finds that the hole density 
modulation is caused by the quantum interference pat¬ 
tern of the phase string effect as a bulk property, which 
cannot be reduced to a conventional standing wave due 
to two counter-pro pag ating Bloch waves under the open 
boundary condition 2*4 Even though the ground state 
at a > a c is concomitant with the momentum split¬ 
ting/Fermi surface (point) reconstruction, the Landau’s 
one-to-one correspondence principle nonetheless breaks 
down here. Indeed, it is no longer meaningful to try 
to identify the hole state |'F< 3)111 in Eq. © with a con¬ 
ventional quasiparticle since an adiabatic connection in 
the Landau’s paradigm is broken down in such a trans¬ 
lational symmetry breaking regime. In particular, it has 
been pointed out that the self-localization of the charge 
is inevitable at a > a c based on the variational form of 
Eq. j5j. 

Finally, the similar symmetry breaking state |4 /g)iii is 
in principle applicable to the t-J ladders with more legs. 
It includes the two-dimensional limit, which is relevant 
to the high-T c problem in the cuprates. The generaliza¬ 
tion of the present VMC for the single hole problem is 
straightforward, although the DMRG convergence gets 
more and more difficult with the increase of the leg num¬ 
ber. A VMC study along this line is currently under¬ 
way. Furthermore, the DMRG calculation has shown 18 
a strong binding between two doped holes in the two-leg 
ladder, implying the ground state Eq. ([2]), which can be 
also studied by VMC in the future. 


ACKNOWLEDGMENTS 

We acknowledge stimulating discussions with R. Q. He, 
H. C. Jiang, D. N . Sheng, C. S. Tian, and J. Zaanen. 
Work was supported by the NBRC (973 Program, Nos. 
2015CB921000 and 2011CBA00108), by Tsinghua Uni¬ 
versity’s ISRP, and in part by Perimeter Institute for 
Theoretical Physics. Research at Perimeter Institute is 
supported by the Government of Canada through Indus¬ 
try Canada and by the Province of Ontario through the 
Ministry of Research and Innovation. 


1 P. W. Anderson, Science 235, 1196 (1987) 

2 P. W. Anderson, P. A. Lee, M. Randeria, T. M. Rice, 
N. Trivedi, and F. C. Zhang, Journal of Physics: Con¬ 
densed Matter 16, R755 (2004), and references therein. 

3 P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006), and references therein. 


4 B. Edegger, V. N. Muthukumar, and C. Gros, Advances 
in Physics 56, 927 (2007), and references therein. 

3 D. N. Sheng, Y. C. Chen, and Z. Y. Weng, Phys. Rev. 
Lett. 77, 5102 (1996) 

0 Z. Y. Weng, D. N. Sheng, Y.-C. Chen, and C. S. Ting, 
Phys. Rev. B 55, 3894 (1997) 















12 


' K. Wu, Z. Y. Weng, and J. Zaanen, Phys. Rev. B 77, 
155102 (2008) 

° L. Zhang and Z.-Y. Weng, Phys. Rev. B 90, 165120 (2014) 
9 Z.-Y. Weng, New Journal of Physics 13, 103039 (2011) 

10 Y. Ma, P. Ye, and Z.-Y. Weng, New Journal of Physics 
16, 083039 (2014) 

11 P. W. Anderson, Phys. Rev. Lett. 64, 1839 (1990) 

12 S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

13 Z. Zhu, H.-C. Jiang, Y. Qi, C. Tian, and Z.-Y. Weng, 
Scientific Reports 3, 2586 (2013) 

14 S. Schmitt-Rink, C. M. Varma, and A. E. Ruckenstein, 
Phys. Rev. Lett. 60, 2793 (1988) 

15 C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 
6880 (1989) 

1D G. Martinez and P. Horsch, Phys. Rev. B 44, 317 (1991) 

17 Z. Liu and E. Manousakis, Phys. Rev. B 44, 2414 (1991) 

18 Z. Zhu, H.-C. Jiang, D. N. Sheng, and Z.-Y. Weng, Scien¬ 
tific Reports 4, 5419 (2014) 

iM Z. Zhu and Z.-Y. Weng, ArXiv e-prints (2014), 
arXiv: 1409.3241 [cond-mat.str-el] 

20 Z. Zhu, C. Tian, H.-C. Jiang, Y. Qi, Z.-Y. Weng, and 
J. Zaanen, Phys. Rev. B 92, 035113 (2015) 

21 S. R. White, D. ,J. Scalapino, and S. A. Kivelson, Phys. 
Rev. Lett. 115, 056401 (2015) 

22 W. Marshall, Proc. R. Soc. Lond. A 232, 48 (1955) 

23 E. Dagotto, J. Riera, and D. Scalapino, Phys. Rev. B 45, 
5744 (1992) 

23 M. Sigrist, T. M. Rice, and F. C. Zhang, Phys. Rev. B 
49, 12058 (1994) 

23 H. Tsunetsugu, M. Troyer, and T. M. Rice, Phys. Rev. B 
49, 16078 (1994) 

“ !D H. Tsunetsugu, M. Troyer, and T. M. Rice, Phys. Rev. B 
51, 16456 (1995) 

2 ' E. Dagotto and T. M. Rice, Science 271, 618 (1996), and 
references therein. 

28 M. Troyer, H. Tsunetsugu, and T. M. Rice, Phys. Rev. B 
53, 251 (1996) 

C. A. Hayward and D. Poilblanc, Phys. Rev. B 53, 11721 
(1996) 

30 T. K. Lee and C. T. Shih, Phys. Rev. B 55, 5983 (1997) 

31 S. R. White and D. J. Scalapino, Phys. Rev. B 55, 6504 
(1997) 

Y. L. Lee, Y. W. Lee, C.-Y. Mou, and Z. Y. Weng, Phys. 
Rev. B 60, 13418 (1999) 

13,3 J. Oitmaa, C. J. Hamer, and W. Zheng, Phys. Rev. B 60, 
16364 (1999) 

31 F. Becca, L. Capriotti, and S. Sorella, Phys. Rev. Lett. 
87, 167005 (2001) 

130 S. Sorella, G. B. Martins, F. Becca, C. Gazza, L. Capriotti, 
A. Parola, and E. Dagotto, Phys. Rev. Lett. 88, 117002 
( 2002 ). 

■ 3 ° S. Liang, B. Doucot, and P. W. Anderson, Phys. Rev. 
Lett. 61, 365 (1988). 

13 ' Here we express the “bosonic” VB states by electron cre¬ 
ation operators. 

38 J. Zaanen and B. J. Overbosch, Philosophical Transactions 
of the Royal Society of London A: Mathematical, Physical 
and Engineering Sciences 369, 1599 (2011) 

Ja Z.-Y. Weng, Frontiers of Physics 6, 370 (2011) 

40 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983) 

41 A. W. Sandvik and H. G. Evertz, Phys. Rev. B 82, 024407 

( 2010 ) 

^ I. M. Lifshits, S. A. Gredeskul, and P. L. A., Introduc¬ 
tion to the theory of disordered systems (Wiley-Interscience 


Publication, NY, 1988). 

43 R. Sun, L. Zhang, and Z.-Y. Weng (unpublished). 

44 Q.-R. Wang and P. Ye, Phys. Rev. B 90, 045106 (2014) 







13 


Appendix A: VMC for single-hole wave function 


To provide the necessary notations and make this paper more self-contained, we first present the VMC procedure 
for the (half-filled) RVB state following Ref. ([36]) and (|4T]l . Whereafter, the VMC formulas for the single-hole wave 
function Tc) are derived. 

The normalization of RVB state Eq. ( flO| is given by 

(RVB | RVB) = 5>„»„{</|»>. (Al) 

v,v' 


Since w v >w v (v\v) is positive, we can interpret it as a distribution function. The average value of a physical quantity 
O is 


(O) 


(RVB|0|RVB) 
(RVB | RVB) 


y- w v >w v (v'\v) (y'\0\v) 
(RVB|RVB) (v'\v) 


(A2) 


The quantity (v'\0\v)/(v'\v) to be averaged in VMC is usually of order one. For O = S* • S j, we have 


('c'IS, • Sj|u) 

(v'\v) 



3 

4’ 


(A3) 


where i5*°° p = 0 or 1 indicates whether the two sites i and j belong to the same loop in the transposition-graph of 
dimer covers v,v'. 

The most time-consuming part of VMC is the loop tracing in calculating the overlap {t/|u). One way to circumvent 
this problem is to sample the overlap in Monte Carlo by introducing an Ising configuration a (we use a instead of {a} 
for simplicity), besides the two dimer covers v and To combine the VB state and the Ising basis, we introduce 
the notation 


|cr)(cr|u) = 5 v ^\v, a) = r] v ,a\cr), 


(A4) 


where 5 v ^ a = \r] Vi<T \ and r] Vi<7 = (cr|v) = 0, ±1 is zero or the Marshall sign for the ground state wave function of 


antiferromagnetic Heisenberg model. Now the VB state |r>) and the RVB state Eq. (10) can be expressed as 

\v) = ’^ 2 , 5 v ,„\v,<j), 

<7 

I RVB) = y^/w v \v) = Y. 5 v ,crl 


T w v v,a 


(A5) 

(A6) 


The summation is constrained in the space where the dimer cover v and Ising bases a are compatible, i.e., <5„ j(T = 1- 
By using the fact 


N 


(u , |u) = 2 
{v',ct'\v,<t) = 6 a ^>, 


(A7) 

(AS) 


the norm of RVB state is now 


(RVB | RVB) = Y, w v' w v( V '\ v ) = Y, 5 V ',crS Vt aW v ’W v . 

v,v' V,v',(J 


(A9) 


Here, N^F is the number of loops in the transposition-graph of dimer covers v,v'. Note that the right hand side 
of Eq. (A7) is exactly the number of Ising bases compatible with both the dimer covers v and v'. As a result , we 
can sample the mutual compatible triad {v,v',a} configuration space to get the expectation value (O) in Eq. (A2) 
without explicitly calculating the overlap (v'\v): 


,A\ _ sr a a Wv ' Wv ( v '\o\ v ) 

[ ^ (RVB|RVB) (v'\v) 

V,V,(T 1 

(E vy,a Sv '^ S v,<y) W v'W V ■ W\0\v)/{v’\v) 

W V 'W V 

The same trick is used in VMC simulations of the single-hole wave function |^g). 


(A10) 
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1. Single-hole wave function 


We introduce the single-hole “VB” states by removing a spin (an up spin without loss of generality) at site h from 
the half-filled VB states: 


I h, v ) = b h t|v) = 55 6 V}Crh | h, v, a h ), 


(All) 


where ah is an Ising basis on the lattice without site h, while the dimer cover v covers the whole lattice. And by 
analogy with Eq. (A41, we use the notation 


\h,a h ){h,a h \h,v) = 5 v ^ h \h,v,a h ) = r] Vt<Th \h,a h ). 


(A12) 


S v>ah is zero whenever v and ah is not compatible, i.e., ah{i) = Vh(j) for some dimer (i,j) £ v , or the spin of the 
site h! originally connecting the hole site h is not a down spin: ah{h') ^\ r . And again, r] Vjah = 0,±1 is zero or the 
Marshall sign. 

The variational single-hole wave function is obtained from the RVB state Eq. (101 by removing a spin, accompanied 
with a unitary transformation A: 


I^g) = Ay^y^(fe)cfa t |RVB) = y 6 V!<Th A(h,a h )ip h (h)w v \h,v,a h ). 

h h,v,ah 


(A13) 


Here, the hole wave function is normalized as JV Wh{i)\ 2 = 1- The U(l) phase factor A(h,ah) is a function of 
the hole position h and spin configuration ah, and is defined by 


A| h, v, a h ) = A (h, a h )\h , v, a h ) = P A (h, l, a h (l)) \h, v, a h ). 

l^h 


(A14) 


The phase factor A (h, ah) is fractionalized to Yli^h ^ a h(l)) in the last step of Eq. (A14), which is similar to the 
fractionalization of amplitude in the Liang-Doucot-Anderson type RVB state®) We have different choices of phase 
factor A (h, l, a): 

(i) If we choose A (h, l, a) = 1, then A = 1, and I’l'c) and |4 >g) are the same. 

(ii) For 2D rotational invariant system, we choose A(h,l,a) = , where cf> h i = Im ln(^ — zi ) and Zj = Xj + iyj 

is the complex coordinate of site j. 

(iii) For ladder system and A defined in the mainbody of the paper, we choose an anisotropic phase factor charac¬ 
terized by A: A (h , l, a ) = e l9hl5 ” i , where and e l6hl are in the same quadrant and tan 0m = A tan (j>hi. 

Similar to Eq. (A9), the normalization of the single-hole wave function is given by 


(^g|^g) = 5151 ( 55<Wh<Wh I \<Ph(h)\ 2 w V 'W v 




N'°°v -1 


W V 'W V 


due to the inner product of our bases 


— J ^ ^ ,cr° J cy'W v 'W v 

\V,V f ,< 7 ° 


jy loop_i 


(h',v'\h,v) = 8 h ,h'2 
{h , v , ah' \h, v, ah) — 8h,h',a' h • 


(A15) 


(A16) 

(A17) 


These equations bear a resemblance to Eq. (A7) and Eq. (A8). Note that there is a —1 in the exponent on the right 
hand side of Eq. (A16), for the transposition-graph loop containing the hole contributes only one Ising configuration 

jyl°°P_i 

rather than two. The number 2 is also exactly the number of Ising bases compatible with hole position h and 

dimer covers v, v'. In the last line of Eq. (A15), we put an additional summation over the Ising bases cr° on the whol e 
lattice (without hole), such that the configuration space {(u, v', er 0 )} is the same as the half filled case in Eq. (A9|. 
Comparing Eq. (A15) to Eq. (A9), we find that the normalization of '•I'c) and |RVB) are related by 

1 


<*g|*g> = ^(RVBIRVB). 


(A18) 
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2. Superexchange energy 


We start with the expectation value of the Heisenberg superexchange terms, which are easier in the sense that they 
do not change the hole position. The average superexchange energy between sites i and j is 


where 


i H ij) = Y H H S v,vJv'y h ' Wh(h)\ 2 • Re (A* (ft, a' h )A(h, a h )) ■ (h,v',a' h \H^\h,v,a h ) 

v,v' ' ' h a h ,a' h 


(j2v,v\a° S v,^S v > t<7 o^ \w V ’W v • E^(v,v') 

(^^/V,V r ,(7° ^v' ,CT 0 ^ 2^ v '^ v 


E ij{v,v') =J2\ip H (h)\ 2 ES j (h,v,v'), 

h 

r , , (h.v' ,cr' h \H/Ah,v,<Jh) 

E ij(h,v,v)= Y 6 v,a h S v >y h • Re (A (h,a h )A(h,a h )) - ^ ^ ^ -■ 


(A19) 

(A20) 

(A21) 


We can now interpret (\/2)w v 'W v / g\^g) as a probability function in the space of compatible configurations 
(v,v',a°). The superexchange energy can be calculated by averaging Ef :j (v. v r ) in standard Monte Carlo proce¬ 
dure. For a fixed configuration, the energy to be averaged E^(v,v') has summations over hole position ft. and spin 
configuratio ns aj, , v' h ■ 

The Eq. (A21) can be further simplified. The element (ft, v', a' h \H^ |ft, v, ah) makes the summations over spin 
configurations easier as it forces the spin configurations 07 , and a' h are almost the same except on sites i and j . Now 
fix the conhguration (v,v',a°) and the hole position ft, we should distinguish three different situations to simplify 
Eij(h,v,v'): 

(i) Hole site ft coincides with sites i or j, then E^(h,v,v') = 0. 

(ii) Sites i and j belong to different loops in the transposition-graph of dimer covers v, v' . For terms S+S~ and S+ 
in H‘lj , the expectation values are always zero because of the compatibility of the dimer covers and spin configurations 
(one closed loop can not have a single antiferromagnetic domain wall). For diagonal term SfS?, however, although the 
expectation over a fixed spin configuration is not zero, the summation of these terms is zero due to the independence 
of a h (i) and a h (j). 

(iii) Sites i,j belong to the same loop in the transposition-graph. If this loop does not contain the hole site, one 
can show Eq. dA2l| ) becomes 


E?Ah,v,v’) = 2 x 2<°'- 2 • Re (A A/) 


= ^Re(AA 


- J/2 + 2 a O’ _1 


N 10 °P_1 


-J/ 4 

jV lo °P-l 


(A22) 


where the hrst term in the first line comes from the equal contributions of S+ S • and S i SJ , and the second term 


comes from Sf S'- . The result Eq. (A22) is still valid, when the sites i,j , ft all belong to the same loop (but these three 
sites are different). But the origins of each term are different: only S+ S~ or S~Sj~ contributes to the first term. 


Now we turn to the definition of phase factor change A A fj{h) in Eq. (A22), which comes from the phase difference 
between the bra and ket of the single-hole wave function: 

(i) When we choose A(ft, l , a) = 1 by ignoring the phase factor A in the single-hole wave function, the phase factor 
change AA.ft (ft) = 1 and Eq. (A22) becomes E^(h,v,v') = —3J/4 when i,j do not coincide with ft and belong to t he 
same loop in the transposition-graph, which recovers the result of the half-filled energy expectation value Eq. (A3). 

(ii) For a generic fractionalized phase factor A(ft, ah) = Ilz^h A(ft, l, &h(l)), the phase factor change is 


AA(ft ) = A*(ft, i, t)A* (ft, j , |) A(ft, i, 4.)A(ft, j, f). 
Rotational symmetry (A (h,l,a) = e l ^ hl5 ‘ 7 ^) simplifies the above result to 


(A23) 


AAf(ft)| = e l ^ hi = exp ( i Imln 


Zh - Zj 
Zh - Zj 


(A24) 
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(iii) For ladder system with A(h,l,<r) = e l6hl ° ai and tan#/,,; = A tan ((>/,/, we can use Eq. (A23) and have 


aa £(/0 


ladder 


_ e i(8hi-8hi)' 


(A25) 


In summary, the total superexchange energy is ev aluat ed in Monte Carlo by calculating Eqs. (A20), where 
Ef ;j (h. v, v') is zero in some conditions or given by Eq. (A22) otherwise. Same as the (half-filled) Heisenberg model, 
the Monte Carlo configuration space is spanned by two VB states and one spin configuration (v,v',a°) which are 
compatible, with non-negative weight w v 'W v j 2 . 


3. Hopping energy 


Now turn to the expectation value of the hopping term, which moves the hole from one site to another. Direct 
calculation shows: 

( H t) = (VE/A'l/g) X! H ' h , ■Re(ipUti)v h (h)A*{ti,a' h ,)A(h,a h )) 

v,v' ' / h,h' ah,o-' h , 

r / / / 


(ti,v',a h ,\H t \h,v,a h ) 

E V,v',a0 $ V ,<T°J V ', C t°) \w v 'W v ■ E t {v,v') 


(E v,v',cr 0 


2 Wy'Wy 


where the averaged quantity in VMC is 

E t (v,v') = ’Y' j E t (h,h',v,v'), 


h,h' 


E t {h, h', v, v') = ^2 & v ,oJ v 'y h , ■ Re{ip* h (ti)<p h (h)A*(ti , a' h ,)A (h, a h )) • V,afl - 


(A26) 

(A27) 

(A28) 


Now we would like to simplify Eq. (A28). The summation over hole positions h and h! in Eq. (A27) has a constraint 
that h and h! must be neighbouring sites, otherwise the hopping energy vanishes. Now fix a Monte Carlo configuration 
(v,v',a°) and sites h,h’, there are two different situations according to whether sites h and h! belong to the same 
loop in the transposition-graph of dimer covers v and v': 



FIG. 9. (color online). Topological graph of hole down-spin exchange process (sites h and h! belong to the same loop). Black 
and red dots represent sites of different sublattices. Dashed line means the length of the loop is arbitrary, (a) Configuration 
before hopping. Dimers belong to vt (v) are illustrated by black (red) bonds, (b) Configuration after hopping. Spins on this 
transposition-graph must be alternating in the new loop, with hole understood as an up spin. 


(i) h, b! belong to the same loop ■ The only possible incident under the action of H t is the exchange of a 
hole and a down spin (see Fig. [9]). The orthogonal property Eq. (A17) has several constraints: the bond (ij) must 
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be the same as (hh'}; the spin configurations ah and a' h , must satisfy relations ah(h') = a' h ,(h) =4,, ah(l ) = cr//(?)> 
for l ^ h,h!\ spin configuration on sites belong to Lh,h' is uniquely determined by (h, h', v , v'), while for every other 
loop L ^ Lh,h' 1 there are two possible spin configurations. For any given pair of initial and final states with nonzero 
hopping energy contribution, we have (h 1 , v', a' h ,\H t \h, v, ah) = —otijt (the Marshall sign difference of the initial and 


final state cancels the fermion permutation sign). Therefore, the total hopping energy Eq. (A28) is given by 


E hi(h, h', v, v') = -aijt • Re ti, v, v')) , 


(A29) 


where the total phase factor change AAL(/i,/i , ,v,i/) is the only nontrivial value to be calculated. Follow Eq. (A14) 
and the spin configuration constraints, the fractionalized phase factors can be divide into three parts: l = h or h , 
l G Lh,h r , l € L ^ Lh,h'- Accordingly, the total phase factor change AA ^(h, h',v,v') in the hopping process is a 
product of three phase factors: 


AA t H (h,ti,v,v')=A*(h',h,l)A(h,h',l)- ( H A*(h',l,a h (l))A(h,l,a h (l)) 

K l€L hih , 

n \ e n^woMM.^o)) ■ 

\L^L hi h' <Jh(L)=±llGL J 


(A30) 


The spin configuration a h (l) for l £ L ht h‘ i n the second phase factor is determined by (h,h',v,v'). While the 
summation y) fT; (l)=± 1 the third phase factor sums over two possible spin configurations in a loop L different from 

jyloop_ 

Lh,h'- Note that this summation totally gives us 2 oX terms as the value of (h,v'\h,v) (cf. Eq. (A16)), which 
becomes the factor 1/2 when we move it through the product operator of different loops in the third phase factor. If 
we choose A (h,l,a) = e l9hlS,,i for ladder system, the phase factor change becomes 

( \ 

AA^(/i, h,', u, n')| ladder = — [] 

\ l^h,h' 

where the minus sign comes from the first phase factor A*(h',h,\.)A(h,h',1), and A and B denote two sublattices. 


n \ n+n 


0 i(dhl-6 h u) 


(A31) 


/ 


LjtL hth , 


leL leL 
\ieA leB/ 





(a) 


(b) 


FIG. 10. (color online). Topological graph of hole up-spin exchange process (sites h and h! belong to different loops), (a) 
Configuration before hopping, (b) Configuration after hopping. Spins on this transposition-graph also must be alternating in 
the new loop, with hole understood as an up spin. 


(ii) h , h? belong to different loops Lh and Lh>- In this case the hole can only exchange with an up spin (see Fig. 101. 
The final results are parallel to the first case: replace spin down by spin up; replace loop Lh t h’ by loop Lh or Lh> ■ 
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The total hopping energy is 


Otiit 


E{\(h,ti ,v,v') = • Re (ipl{ti)v h (h)&k\^(h,ti,v,vf)) . 


(A32) 


Note that there is fermion permutation sign but no Marshall sign difference in the up-spin hopping process. Therefore, 
there is a minus sign in front of Eq. (A32) comparing to Eq. (A29l. The phase factor change for general A (h,l,a) 
and for the ladder system are 


AA t M (h,ti,v,v') = A*(fc , ,Ji,t)A(M , ,t)- ( II A*(ti,l,a h (l))A(h,l,a h {l)) 

ilELh or L h f 

1 


II o E HA*(h',l,a h (l))A(h,l,a h (l)) 


\L^Lh,L h i a h (L)—±lleL 


(A33) 


AA l t (h, h\v,v') 


ladder 


n 


o i ( S hl-S h / I )S lThWi 


IdzLh or Lfr 


n * n+n 




7 


L?L h ,L h 


leL leL 
<ieA leB/ 


(A34) 


The minus sign in front of the phase factor change Eq. (A311 disappears in Eq. (A34) because A *(h', h, f)A(/i, h' , t) = 
—A*(h',h,l)A(h,h',l.) = 1 for ladder system. The additional numerical factor 1/2 in the hopping energy Eq. (A32) 

rloop 


comes from the fact that the overlap of dimer states is (h,v'\h,v) = 2 N ^y 1 , while there are only A^°°, p — 2 loops 
(those different from Lh and L^) whose spin configurations are not determined. 


In summary, the total hopping energy is evaluated in Monte Carlo by calculating Eqs. (A27), (A29) and (A32l. The 
Monte Carlo configuration space is also spanned by compatible triad (v, v' , <r°), which is the same as the configuration 
space for (half-filled) pure spin model. 


4. Momentum distribution 


To calculate the momentum distribution n k for ITg), we should consider long range hopping process: 


(Aker) — (c^Cka) — pj ^ 


- 1 


( C io C 3°) ■ 


(A35) 


We denote the average value of long range hopping process by T) Jcr = {c. i(J c la ). Different from the calculation of the 
hopping energy which involves only neighbouring sites, i and j in T 1 , J(J can be the same site or separate far from each 
other. Direct calculation shows: 


Tijo = (cLcj g ) = Y H 5 v,aJv',o' h , ■ Re {ip* h (ti)y h ( h)A*(h', <j' h ,) A (h, a h )) 

v,v' ' ' h,h' ah,cr' h / 

■ (h! ,v' ,<j’ h ,\c\ a c ja \h,v,<j h ) 

(T,v,v',<r° s v,< 7 °jv',a° h ) \w V 'W V ■ T ijrJ (v,v') 


(^2v,v',cr 0 ^v,a°^v'2^ v '^ v 

where Tij a (v,v') is to be calculated in every Monte Carlo measurement step. 

The expressions of Tij a (v, v') are different for i = j and i ^ j. Let us consider i = j first: 

Tiia(v,V ) = y ' ..loop ^ ' &v,a h Sy' ,cr h \Vhjti) \ 
h^i ^ ’ (Th 

= y Wh(ti )\ 2 kr latt + -m 1 - C latt ))+(! - 

h^i ^ 


(A36) 


(A37) 
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Here, d'° } op = 0 (1) when sites i and j belong to the same loop (different loops) in the transposition-graph. Similarly, 
<5|“ blatt denote whether sites i and j belong to the same sublattice. In fact, T iia is the occupation number n tcr , which 
can be used to calculate the hole density n/: 


n< 


= 1 - ^n io = ! - ^2{T iia (v,v')) v y = 1 - ^\Vh{h)\ 2 ^°° p + (1 - <^°° p ) = \v>h{i)\ 


h^i 


(A38) 


If we use the normalization JA \ip h (i)\ 2 = 2 as in the main body of the paper, then the above result is We 

conclude the hole density is simply \iph(i)\ 2 , which can be easily calculated without VMC. 

On the other hand, for i j, we have 


T lja {v,v') = Afloo ^ 1 5 v,<rjv',*> hl ■ Vh(ti)v h (hW(ti,a' h ,)A{h,a h ) ■ Vcy h V^ h 


pjioop 
O iV V' 

z ’ h,h' <7 h y h , 


■ <0| c N(T , n ■ ■ ■ Ci ffJ c] l , t ( 4 ff c j(T )c /lt c p CT , • • • c^JO) 

1 


J J 


Ar ioop_ 
2' 


inY ■ w* h {j)vh(iW{j,a'j)A{i,(Ti) ■ rjfj i rj a 




(A39) 


The minus sign in the last line is the fermion sign which comes from the permutation of the c and c 1 operators. There 
are 2 3 = 8 different cases in which the final result of Elj a (v, v') takes different forms. They are classified according to: 
(1) the spin er; (2) whether sites i and j belong to same sublattices; (3) whether sites i and j belong to the same loop 
in the transposition-graph of dimer covers v and v'. These results are summarized in Table [IJ Note that the signs in 
front of the results are combinations of the fermion signs and the Marshall signs. The factor 1/2 stems from the fact 
sites i and j belong to different loops. 


TABLE I. Summary of Tij a (v, v') (i y£ j) in 2 3 


8 different cases. The phase differences AAL are given by Eqs. (A401-( A451. 


cases 

a 

sublattices 

loops 

Tij a (v,v') 

i 

t 

different 

same 

0 

2 

t 

different 

different 

-^hU)<Ph(i) AAj£ 

3 

t 

same 

same 

-<PhU)<Ph(i) AA f/ 3 

4 

t 

same 

different 

-l‘PhU)‘Ph( i ) AA ija 

5 

4 

different 

same 

< PhU)‘Ph{i) AA 

6 

4 

different 

different 

0 

7 

4 

same 

same 


8 

4 

same 

different 

0 


All the phase difference AA fj a in the last column of Table |l] can be divided into three parts: 

AA = AAi ■ AA 2 • AA 3 . (A40) 

The expressions for the above three phase difference parts are: 

(i) AAi comes from terms l = i or 1= j: 

AAi = A*(j,i,a)A(i,j,a) = 6^ t - d^. (A41) 

(ii) AA 2 comes from terms where site l (/ y^ i,j ) belongs to the same loop as site i or j in the transposition-graph 
of v and v'\ 


AA 2 = H A*(j,/,a i (Z))A(t,Z, f 7 i (0). 

l(zLi or Lj, 


(A42) 
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Note that the spin configuration on site l (l £ Li or Lj) is totally fixed in each of the eight cases (for instance, see 
Fig. [ 9 ] for case-5, and Fig. [To] for case-2). If we choose A (h,l,<j) = e l(lhlSa> ' for the ladder system, the phase difference 
becomes 


AA' 


21 ladder 


n 




(A43) 


l(zLi or Lj, l^i,j 


(iii) AA 3 comes from terms where site l belongs to different loops as site i or j in the transposition-graph of v and 
v'. The spin configurations on these loops have two possibilities. 

aa 3= n \ e n A *(?d^i(o) A (id,^(o)- ( A44 ) 

L^Li ,Lj cn (Z/)=i l£L 

Similarly, for the ladder system, this phase difference is given by 


AA. 


3 I ladder 



n 


i(8a-6 


(A45) 


l£L,l£B 1 


The spin configuration on each loop L (L ^ Li,Lj) has two possibilities ( cri(L ) = ±). AA 3 on each loop is obtained 
by averaging the phase differences of these two possibilities. 

The hopping energy calculation in Sec. |A 3| can be viewed as special cases in Table[IJ The up-spin hopping on nearest 
bond corresponds to cases 1 and 2 in this table. Since case-1 has zero result, only case-2 contributes to TL^.. This is 


exactly Eq. (A32) if —aijt is added. Similarly, for down-spin hopping, case-5 gives the result Eq. (A29), while case-6 
has no contribution. 

To sum up, the momentum distribution n k is calculated from Tij a Eq. ( A36|. The eight cases of Tij a (v,v') are 


summarized in Table|lJ with the phase differences AA Jj a given by Eqs. (A40|-(A45). 


5. Quasiparticle weight 

J I VTr —. \ 12 


The quasiparticle Z k is defined by |(RVB|cL.|’J , g )\ with normalized I’I'g) and |RVB), or equivalently 


^k = 


|(RVB|cj; t |vI/ G )| 2 

(RVB|RVB)(4' G |4' G ) 


= 2 


Vs* 


-ik-R, 


(RVB|cL| 4' g ) 


(RVB|RVB) 


= 2 


-L y 


e- lkR 'z t 


(A46) 


where the normalization relation Eq. (A18) is used. Zi is roughly the average of phase string factor e ,n * and defined 
by 


where 


(RVB|cie- iA -c lT |RVB) 


Zi = 


A 


(RVB|RVB) 

(j2vy,cr0 S v,<r°8v',< t °) W v >W v ■ Zi(v,v') 


,cr 


0 ^V } (T a ^v\cr^ j Wy'Wy 


(A47) 

(A48) 


= ip h {i) ( i ( n \ n A (bU(0)]. 

Z a(L)=±lGL 


(A49) 


Similar to the VMC simulation for |RVB), Z; t is ob taine d by averaging Zi(v, v') with respect to (v, v', a 0 ) with weight 
w v rw v . Z k is then calculated directly by using Eq. (A46). 
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Appendix B: Sign structure of the a-t-J model 

In this appendix, we will show explicitly that the sign structure of the a-t-J model, on a bipartite lattice in arbitrary 
dimension and hole concentration, is the Marshall sigrf 22 !, instead of the phase string for the t-J model. In particular, 
the Bloch-like wave function |ko)eL with ko = (0,0) satisfies the sign structure requirement. 

Let us start with a generic single-hole-doped wave function which is denoted as 



(Bl) 

bW 


Here the basis state is defined as 


IbM) = cqJM), 

(B2) 

where the half-filled Marshall basis is given by 


im> = (- l ) Jvf c t lCT1 cL 2 ---cL N l°), 

(B3) 

with NP the number of down spins belonging to sublattice B. 

The sign structure is determined by the off-diagonal elements of the a-t-J Hamiltonian, 
off-diagonal elements of the hopping terms of the a-t-J model in the basis Eq. (|B2|) are 

Specifically, The nonzero 

~t ( 3, W}\4 t c jt \i, {o-}) = -t ({o-'}|c^c] t c it c 4 |{er}) 


= t(K}|S+57|M)<0, 

(B4) 

t {j, W) lctj_Cj4.fi, M) = t ({<T , }|ct J _c] J _Cj4.c i 4,|{<r}) 


= -t ({a'jlmnj |{cr}) < 0. 

(B5) 

On the other hand, the nonzero off-diagonal elements of the superexchange terms are 


J - (/l,{cr'}|S'+5'"|/l,{cr}) < 0. 

(B6) 


The nonnegativity of both Eq. (B4) and Eq. (B6| is owing to the Marshall sign (—1)^ in the basis Eqs. (B2) and 
|. We conclude the off-diagonal elements of the a-t-J Hamiltonian are all non-positive in the basis Eq. (B2|. 


As a result, according to the Perron-Frobenius theorem, the ground state of the a-t-J model has the form of Eq. (B1) 
with ip(i, {er}) > 0. That means, the sign structure of the a-t-J model is exactly the Marshall sigrP^ in the basis 
Eq. (B3|, the same as the Heisenberg spin model. In particular, if we ignore the spin polaron effect, the Bloch-like 


wave function |k 0 )BL at k 0 = (0,0) should well describe the a-t-J model, for it satisfies the sign structure of this 
model. Indeed, Fig. [8] illustrates the overall agreement of the Bloch state |k 0 )BL at k 0 = (0, 0) with the DMRG result. 













